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Abstract: We have proposed that the abiogenesis of life around the beginning of the Archean may 
have been an example of “spontaneous” microscopic dissipative structuring of UV-C pigments under 
the prevailing surface ultraviolet solar spectrum. The thermodynamic function of these Archean 
pigments (the “fundamental molecules of life”), as for the visible pigments of today, was to dissipate 
the incident solar light into heat. We have previously described the non-equilibrium thermodynamics 
and the photochemical mechanisms which may have been involved in the dissipative structuring of 
the purines adenine and hypoxanthine from the common precursor molecules of hydrogen cyanide 
and water under this UV light. In this article, we extend our analysis to include the production of 
the other two important purines, guanine and xanthine. The photochemical reactions are presumed 
to occur within a fatty acid vesicle floating on a hot (~80 °C) neutral pH ocean surface exposed to 
the prevailing UV-C light. Reaction—diffusion equations are resolved under different environmental 
conditions. Significant amounts of adenine (~10-5 M) and guanine (~10-© M) are obtained within 
60 Archean days, starting from realistic concentrations of the precursors hydrogen cyanide and 
cyanogen (~1075 M). 


Keywords: origin of life; dissipative structuring; non-equilibrium thermodynamics; prebiotic chem- 
istry; abiogenesis; adenine; guanine; hypoxanthine; xanthine; purines; hydrogen cyanide; cyanogen 


PACS: 05.70.Ln; 91.62.Uv; 91.62.+g; 91.62.Kt; 91.62.Np; 91.70.h; 87.14.¢k; 87.14.¢n; 87.15.-v; 87.16.Dg; 
87.23.-n 


MSC: 92C05; 92C15; 92C40; 92C45; 92D15; 80A99; 82C99 


1. Introduction 


Possible chemical or photochemical routes to the production of the purines, including 
adenine, hypoxanthine, guanine, xanthine, and 2,6-diaminopurine from common precursor 
molecules have been extensively studied since the first experiments of Miller and Urey 
in 1953. Sources of free energy for synthesis have included heat gradients, lightning, UV 
light, and shock waves. Plausible precursor molecules are hydrogen cyanide (HCN) and 
water, along with the photochemical and hydrolysis products of HCN, which include 
formamide (H2N-CHO), ammonium formate (NH4HCOz2), cyanogen (NCCN) and cyanate 
(OCN  ) [1-13]. These experiments, however, were usually carried out at alkaline pH and 
unrealistically large HCN concentrations (>0.1 M) because, only under these conditions, 
polymerization of HCN to cis-DAMN (an intermediate on route to adenine and guanine) 
prevails over hydrolysis. 

From the perspective of the “Thermodynamic Dissipation Theory of the Origin of 
Life” [14-28], descriptions of molecular synthesis alone, however, are insufficient to capture 
life’s vitality, which includes proliferation, dynamics, and evolution. The vitality of life 
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requires instead a constant dissipation of an imposed thermodynamic potential. Of the 
above listed free energy sources employed in synthesis experiments, only UV light would 
have been sufficiently energy dense and continuously available throughout the Archean. 
Our theory therefore affirms that the “fundamental molecules of life” (those in the three 
domains of life) were, at their origin, microscopic photon dissipative structures (i.e., pig- 
ments) which arose “spontaneously” under the thermodynamic imperative to dissipate the 
prevailing UV-C (210-285 nm) solar photon spectrum (Figure 1) [14,15,21,26]. 
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Figure 1. The spectrum of UV-C light available at Earth’s surface before the origin of life at approxi- 
mately 3.9 Ga and until at least 2.9 Ga (curves black and red respectively) and perhaps throughout 
the entire Archean until 2.5 Ga. CO2 and probably some H2S were responsible for absorption at wave- 
lengths shorter than ~205 nm, and atmospheric aldehydes (common photochemical products of CO 
and water) absorbed between about 285 and 310 nm [29], approximately corresponding to the UV-B 
region. Around 2.2 Ga (green curve), UV-C light at Earth’s surface was completely extinguished by 
oxygen and ozone resulting from organisms performing oxygenic photosynthesis. The yellow curve 
corresponds to the present surface spectrum. Energy fluxes are for the sun at the zenith. The names 
of the fundamental molecules of life are plotted at their wavelengths of maximum absorption; nucleic 
acids (black), amino acids (green), fatty acids (violet), sugars (brown), vitamins, co-enzymes and 
cofactors (blue), and pigments (red) (the font size of the letter roughly correspond to the relative 
size of their molar extinction coefficient). In this work, 125 W m~? um7! is the “nominal” incident 
intensity at 260 nm, but we also consider lesser and greater intensities due to the uncertainties in the 
COy atmospheric concentration at the origin of life [30]. Adapted from Michaelian and Simeonov [18]. 


In this paper we model the origin of life as a photochemical dissipative structuring 
process with reaction—diffusion occurring within a fatty acid vesicle floating at the ocean 
surface, permeable to the precursor molecules HCN, cyanogen and HO, and under a 
UV-C photon spectrum arriving at the Earth’s surface during the Archean (Figure 1). The 
concentrations of the relevant product molecules are determined as a function of ocean 
surface temperature, and as a function of incident light intensity, varying the amplitude 
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but not the shape of the spectrum of Figure 1. Starting from a realistic environmental 
HCN concentration of 6 x 10~° M, we show that adenine and guanine concentrations of 
2.4 x 10-° M and 7.6 x 10-7 M respectively, could have been obtained in only 60 Archean 
days at neutral pH, 80 °C, and under a nominal UV-C flux of 125 W m2 um! at 260 nm. 
These final purine concentrations increase exponentially with temperature (roughly an 
order of magnitude for each 10 °C increase in surface temperature) and have a linear but 
complex dependence on light intensity. 


2. Thermodynamics of Dissipative Structuring 


We have identified the long wavelength part of the UV-C region (~210-285 nm), plus 
part of the UV-B and UV-A regions (~310-360 nm), of the solar spectrum as the thermo- 
dynamic potential which could have driven the molecular dissipative structuring [21,26], 
proliferation, and evolution relevant to the origin of life. This light prevailed at Earth’s 
surface from the Hadean, before the origin of life (somewhere near the beginning of the 
Archean ~3.9 Ga), and for perhaps as long as 1400 million years until the formation of an 
ozone layer at about 2.5 Ga [29,31,32] after natural oxygen sinks (e.g., volcanic reducing 
gases such as CHy, CO, HS and Fe*”) became overwhelmed by organisms performing 
oxygenic photosynthesis, possibly fertilized by volcanic phosphorus [33]. 

The relevance of this particular region of the solar spectrum, corresponding to the 
Archean atmospheric window of transparency, to the dissipative structuring of the funda- 
mental molecules is that longer wavelengths do not contain sufficient free energy to directly 
break double covalent carbon bonds, while shorter wavelengths contain enough free energy 
to destroy carbon based molecules through successive ionization or fragmentation. 

A number of empirical evidences support our dissipative structuring conjecture for the 
fundamental molecules of life. First, the wavelength of maximum absorption of many of 
these molecules coincide with the predicted window in the Archean atmosphere (Figure 1). 
Secondly, many of the fundamental molecules of life are endowed with peaked conical 
intersections [26] giving them broad band absorption and high quantum yield for internal 
conversion, i.e., extremely rapid (picosecond) dissipation of the photon-induced electronic 
excitation energy into vibrational energy of molecular atomic coordinates, and finally into 
the surrounding water solvent [21,34]. The most convincing evidence of all, however, is 
that many photochemical routes to the synthesis of nucleic acids [35], amino acids [36], 
fatty acids [23], sugars [37], and other pigments [18] from common and simple precursor 
molecules have been identified at these wavelengths and the rate of photon dissipation 
within the Archean window generally increases after each incremental step on route 
to synthesis, a behavior strongly suggestive of dissipative structuring in the non-linear 
non-equilibrium regime [21,23,26]. Finally, even small changes (e.g., tautomerizations) or 
additions to the fundamental molecules destroys these special optical properties. 

The photochemistry involved in the non-equilibrium thermodynamics of microscopic 
dissipative structuring and proliferation of the fundamental molecules, and how this 
leads to the evolution of molecular concentration profiles with ever greater solar photon 
dissipation capacity (corresponding to increases in global entropy production) has been 
described in detail elsewhere (Michaelian [26]). Here we emphasize only those aspects 
most important to the dissipative structuring of the purines guanine and xanthine. 


3. The Dissipative Structuring of the Purines 
3.1. The Model 


Cyanogen, HCN, and its hydrolysis product formamide, were first recognized by Ed- 
uard Pfliiger in 1875 to be the probable precursors of life [38]. The fundamental molecules; 
nucleic acids [6,35], amino acids [39,40], fatty acids [37], and even simple sugars [41,42] 
have all been synthesized from these precursors. 

HCN is found throughout the cosmos [43,44]. Astronomical observational data in- 
dicates that HCN and water can be found together in unexpectedly high abundance in 
many hot star forming regions [45]. It may be, therefore, that HCN was incorporated into 


Entropy 2022, 24, 1027 


4 of 26 


the bodies that formed the Earth and was gradually out-gassed to the surface, as is now 
suspected for Earth’s water. Cyanide has, in fact, been found in abundance in carbonaceous 
chondrite meteorites [46]. Furthermore, if Earth’s atmosphere during the Hadean and 
Archean had been weakly reducing [30] (see below), HCN production on Earth could 
have resulted from CH, photolysis by the solar Lyman alpha line (121.6 nm) in the upper 
atmosphere, leading to CH* radicals attacking No, or, UV (145 nm) photolysis of Nz and 
the resulting radicals attacking CH or CH? [47,48]. 

Cyanogen, NCCN or (CN)o, can be generated from HCN either photochemically [49] 
or thermally [50]. The Lyman-a line on atmospheric HCN will produce, through photolysis, 
the CN™ radical with high quantum yield (reaction R177 from Table 1 of Zahnle [51]) and 
this radical gives cyanogen after interaction with a second HCN molecule (reaction R199 
from Table 1 of Zahnle). Both HCN and cyanogen have been found in the atmosphere of 
Titan [48] where they probably derive from similar UV photochemical processes in the 
upper or middle atmosphere. 

In Michaelian [26] we have addressed the issue of the expected low concentration of 
HCN in the Archean oceans leading to a low polymerization rate of HCN with respect to 
its rate of hydrolysis. Traditional approaches to the origin of life have “resolved” this issue 
by invoking eutectic freezing to increase the solute HCN concentration to values sufficient 
for biasing the reactions towards polymerization [52-54]. This low temperature approach, 
however, reduces significantly all thermal reaction rates. Instead, we have emphasized the 
existence of an ocean surface microlayer of high organic concentration, 104 or 10° times 
greater than that of bulk water (see Michaelian [26] for discussion and references). We also 
assume the existence of fatty acid vesicles of ~100 um diameter which would allow the 
incident UV-C light, as well as the small precursor molecules, HCN, cyanogen, and H2O, 
to permeate the vesicle bi-layer wall (Figure 2), while trapping inside the reaction products 
due to their larger sizes and larger dipole moments (Table 1). This allows the intermediate 
and product molecules, as well as the heat from UV-C photon dissipation, to accumulate 
within the vesicle. 


UVC Air 


Figure 2. Fatty acid vesicle floating at the ocean surface microlayer of high organic concentration, 
transparent to UV-C light and permeable to the precursor molecules having small size and dipole 
moments; H2O, HCN (H), cyanogen (Cg), and formimidic acid (Fa), but impermeable to the interme- 
diate and final photochemical products (e.g., ammonium formate (Af), cis- DAMN, AICN, adenine 
(A), guanine (G)), which are larger in size and have larger dipole moments (Table 1). 
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Table 1. Nomenclature, chemical formula, abbreviation (used in the text and in kinetic equations), 
label in Figure 3, wavelength of maximum absorption Ajax (Figure 1), molar extinction coefficient 
at that wavelength €y,,, electric dipole moment 1, and the topological polar surface area (TPSA), 

uae 


of the molecules involved in the photochemical synthesis of the purines. Values marked with 
are estimates, obtained by comparing to similar molecules, for cases where no published data have 


been found. 
i A “A TPSA 
Name scilarg Abbrev. in Text ee Figure3 = Amax nm Preis 4 p#IDI [A2] 
hydrogen cyanide HCN HCN H 1 2.98 23.8 
cyanogen NCCN NCCN Cg 10 0.00 47.6 
formamide H2N-CHO formamide : F 220 60 [55,56] 4.27 [57] 43.1 
formimidic acid = H(OH)C=NH_ formic acid Fa 220 60 114[57] 43.1" 
ammonium formate NH,HCO) ammonium formate Af +/—,2.0* 41.1 
diaminomaleonitrile CyHyNg cis-DAMN (DAMN) Cc 2 298 14,000 [58] 6.80 [59] 99.6 
diaminofumaronitrile CyHyNg ae al if 3 313 8500 [58] 1.49 [59] 99.6 
2-amino-3- 
iminoacrylimidoyl C4HyNg AIAC J 4 275 9000 [5,60] 1.49 99.6* 
cyanide 
4-aminoimidazole-5- ] 
carbonitilé C4HyNg AICN I Z 250 10,700 [58] 3.67 78.5 
4-aminoimidazole-5- * 
carboxamide C4HeN,O AICA L 10 266 [39] 10,700 * 3.67 97.8 
5-(N’-formamidiny])- 
1H-imidazole-4- C5H5Ns5 amidine Am 250 10,700 [61] 6.83 * 80.5 * 
carbonitrileamidine 
adenine CsHs5Ns5 adenine A 8 260 15,040 [62] 6.83 [63] 80.5 
hypoxanthine C5HyN4gO hypoxanthine Hy 9 250 12,500 [64] 3.16 70.1 
xanthine C5H4N4O2 xanthine Xa 12 268 [4] 9300 [2] 4.46 [6] 86.9 [11] 
guanine CsHsN5O0 guanine G 11 252 [62] 14,090 [62] 5.45 [65] 96.2 [65] 


The existence of amphipathic fatty acid hydrocarbon chains, which through Gibb’s free 
energy minimization spontaneously form lipid vesicles at the ocean surface, is a common 
assumption in origin of life scenarios [66,67]. Deamer and Georgiou [68,69] have discussed 
in detail the many roles of lipid membranes in life’s origin. In order to maintain vesicle 
integrity at the high surface temperatures considered here at the beginning of the Archean, 
~80 °C [70-75], these fatty acids would necessarily have been long (~18 C atoms) and 
cross linked through UV-C light, which improves stability at high temperatures and over a 
wider range of pH values [23,76], and provides resistance to salt flocculation. In the anoxic 
Archean ocean surface environment, UV-C light would promote cross-linking between 
consecutive fatty acids rather than cause oxidative damage. 

Archean fatty acids may have been aromatic, giving them improved stability at 
higher temperature for shorter lengths due to a van der Waals stacking interaction be- 
tween the rings of the interior and exterior layers of the bilayer vesicle. There is, in fact, 
a predominance of 16 and 18 carbon atom fatty acids in the whole Precambrian fossil 
record [77,78]. Today, the aromatic fatty acid 10-phenyltridecanoic acid (16 carbon atoms) 
and other w-pheny] fatty acids, such as 12-phenyldodecanoic acid (18 carbon atoms) and 
14-phenyltetradecanoic acid (20 carbon atoms), are found in the highly salt resistant ma- 
rine bacteria V. alginolyticus [79], and w-cyclohexyl] fatty acids are found today in some 
acidophilic-thermophilic bacteria isolated from hot springs [80]. 

From the perspective of our dissipative structuring scenario, aromatic fatty acids are 
also interesting because the conjugated aromatic rings have wide band UV absorption cross 
section within the 210-285 nm region and a conical intersection for rapid internal conversion 
to the ground state, giving them large photon dissipation efficacy and photochemical 
stability under this UV-C light. Indeed, we have suggested that fatty acids, just like the 
nucleic acids and aromatic amino acids, could have been dissipatively structured from 
HCN and CO} saturated water under this UV-C light at moderate temperatures on the 
ocean surface [23] because these optical characteristics are the hallmarks of dissipative 
structuring [21,26]. 
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3.2. The Photochemical Route to the Purines 


Our proposed photochemical route to the purines from the precursors HzO, HCN and 
cyanogen is based on the experimental data of Ferris and Orgel [35] and Sanchez et al. [6], 
as well as the quantum mechanical calculations of Boulanger et al. [60], and is given in 
Figure 3. 


H 
i ‘i » 
H,0 
N 


(9) Hypoxanthine (8) Adenine NHe2 
Amax 250 nm Amax 260 nm _ HCN 
N 
S 
4HCN C _hv_, —_hv, _hv_, L- 
298 nm nm 313 am nm 275 nm nm J 
a 
(1) (2) cis-DAMN (3) trans-DAMN (4) AIAC (5) (6) (7) AICN fi H2,0 
Amax 250 nm 
0. N Sey HN 
YI [> 2 + HCN NCCN r» 
HN 
HN 
(12) Xanthine fe} a : ) Guanine (10) AICA 
Amax 268 nm Amax 252 nm Amax 266 nm 


Figure 3. The photochemical synthesis of adenine and guanine from 4 molecules of hydrogen cyanide 
(HCN) in water, plus a final HCN for adenine (Ferris and Orgel (1966) [35,60]), and a cyanogen 
molecule for guanine (Sanchez et al. (1968) [6]). Four molecules of HCN are transformed into the 
smallest stable oligomer (tetramer) of HCN, known as cis-2,3-diaminomaleonitrile (cis- DAMN) (2), 
which, under a UV-C photon flux isomerizes into trans-DAMN (3) (diaminofumaronitrile, DAFN) 
which by absorbing two more UV-C photons transforms into an imidazole intermediate, 4-amino- 
1H-imidazole-5-carbonitrile (AICN) (7). Hot ground state thermal reactions with another HCN 
molecule, or its hydrolysis product formamide (or ammonium formate), leads to the purine adenine 
(8). Hydrolysis of adenine (8) leads to hypoxanthine (9). Hydrolysis of AICN (7) leads to AICA (10) 
and a hot ground state thermal reaction with a cyanogen molecule leads to guanine (11). Hydrolysis 
of guanine (11) leads to xanthine (12). This is a microscopic dissipative structuring process under 
UV-C light [21,26] which ends in the purines; adenine, hypoxanthine, guanine and xanthine, all 
UV-C pigments with large molar extinction coefficients, at 260, 250, 252 and 268 nm respectively, 
and peaked conical intersections providing rapid dissipation of photons around the wavelength of 
maximum intensity of the Archean surface solar UV-C spectrum (Figure 1). 


The route to guanine results from the hydrolysis of 4-aminoimidazole-5-carbonitrile 
(AICN) to form 2-amino-3-iminoacrylimidoyl cyanide (AICA) (step (7) to (10) of Figure 3) 
which then combines through a thermal chemical reaction with cyanogen (NCCN) to form 
guanine plus a HCN molecule [6]. The reactions on route to adenine and guanine are 
in competition with hydrolysis and UV lysis, and these relative rates are dependent on 
concentrations, temperature, pH, the presence of metal ion- and organic intermediate prod- 
uct catalysts, as well as the wavelength dependent intensity of the incident UV spectrum 
(Section 4). Some of these complexities have been studied by Sanchez et al. [5,6]. 

In the following subsection we present our simplified out-of-equilibrium kinetic 
model for 5HCN — adenine and 4HCN + H2O + NCCN — guanine + HCN (Figure 3) 
photochemical reaction system occurring within a fatty acid vesicle floating within the 
high organic concentration microlayer of a hot (~80 °C [70—-75]) Archean ocean surface 
under the UV surface spectrum of Figure 1. Favorable comparisons of our simulated 
photochemical production rates to experimental data have been given in Michaelian [26] 
(this Special Issue). 
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3.3. The Kinetic Equations 


The kinetic equations for the chemical and photochemical reactions are resolved 
numerically. For this perturbed non-linear reaction—diffusion system it will be shown 
in Section 4 that various stationary state solutions exist with spatial symmetry breaking 
resulting in the highest concentration of purines at the center of the vesicle. This would 
then facilitate a subsequent UV-C polymerization of nucleobases into short oligos (of 
course assuming the possibility of UV-C-assisted synthesis of ribose from similar precursor 
molecules [42,81], a UV-C induced nucleoside formation from a base and ribose [82], and a 
temperature [83] or formamide-catalyzed [84] phosphorylation of the nucleosides [85,86], 
but these reactions will not be considered here). This stationary state coupling of reactions 
to diffusion, leading to symmetry breaking (i.e., particular regions of high concentration of 
the products), was shown to occur for purely thermal reactions with different activator and 
inhibitor diffusion rates by Turing [87] and studied more generally as dissipative structures 
under the framework of Classical Irreversible Thermodynamic (CIT) theory by Glansdorff 
and Prigogine [88]. 

Nomenclature, chemical formula, and abbreviations used throughout the text for 
the concentrations of the participating chemical species involved in the photochemical 
reactions leading to the purines shown in Figure 3 are given in Table 1. Their photon 
extinction coefficients, polar surface area and dipole moments (related to permeability of 
the vesicle wall), are also given. 

From an analysis of published experiment [1-13,35,55,56,58,89-95] and time-dependent 
density functional calculations [60,96], the chemical and photochemical reactions listed 
in Table 2 will occur in the photochemical dissipative structuring of the purines from 
HCN, H20, and NCCN. Below, we describe in detail only those reactions pertinent to the 
production of guanine and xanthine from AICN. Reactions leading to AICN and adenine 
have been described in detail elsewhere (Michaelian [26]). 


Table 2. Reactions involved in the photochemical synthesis of the purines (see Figure 3). Temperature 
T is in Kelvin and all kinetic parameters were obtained at pH 7.0. 


# Reaction Reaction Constants 

1 H ky F ky = exp(—14,039.0/T + 24.732); gol hydrolysis of HCN [5,91,94] 

2 Y220+ F — Fa q2 = 0.05 [55,56,90,92,93] 

3 ‘229+ Fa + H+ H2O 93 = 0.03 [92,93,97] 

4 F ky Af kg = exp(—13,587.0/T + 23.735); s~!; hydrolysis of formamide [93,94] 

5 4H oc ks = 1/(exp(—AE/RT) + 1) - exp(—10, 822.37/T + 19.049); M-! s!; AE = 0.61 kcal mol™! [5] 
6 an 7 ke = 1/(exp(+AE/RT) +1) - exp(—10,822.37/T + 19.049); M—! s—1; tetramization[5] 
7 4H Tes 20a Tacs kz = (1.0/(1.0- 0.01)) exp(—(10, 822.37 — 728.45) /T + 19.049); M~? s~! [5] 

8 4H +T+Cg “S2T+Cg kg = k7;M~? s“! [5] 

9a 298 + CoT 9g = 0.045 [58] 

9 = 9313 +T 9 C or = 0.020 [5,58,60] 

10 313 + T > J io = 0.006 [5,58,60] 

11 yo75 tJ I qq = 0.583; T > I; fo X 4 = 0.0034 [5,60] 

1200 pkey ky = exp(—E,/RT + 12.974); s~1; Eg = 19.93 kcal mol~!; hydrolysis of AICN [6] 

13 oppy ag AGRE kig = exp(—E,/RT + 12.973); M-! s~1; E, = 6.68 kcal mol! [98,99] 

14 LE + Fa “4 Am + Fa +H,O ky = exp(—E,/RT + 12.613); M~! s~!; E, = 19.90 kcal mol! [100] 

15 7250 + Am>A 415 = 0.060 [61] 

16 A kg Hy kg = 10(—5902/T+8.15). 5-1. valid for pH within 5 to 8; hydrolysis of adenine [89,95] 

17 L+Ceg Neo4H ky7 = exp(—E,/RT +15.52); M~'s—!; AICA + Cyanogen, E, = 18.49 kcal mol! [6] 
18 G kig Xa kig = 10{~6930/T+9.40) . 5-1. valid for pH within 5 to 8; hydrolysis of guanine [89,95] 
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Table 2. Cont. 


# Reaction Reaction Constants 
19 298 + CoC qi9 = 0.955 
20 -433+ToT 420 = 0.972 
21 7275 + J- J qu = 0.417 
22 250 + Am + Am 922 = 0.940 
23 yo50 t1 I 423 = 1.000 
24 7266 + LoL 924 = 1.000 
25 260 + AOA 925 = 1.000 
26 7250 + Hy > Hy 926 = 1.000 
27 252 + GG 927 = 1.000 
28 26g + Xa —> Xa 928 = 1.000 


The following is a description of the reactions given in Table 2 by reaction number. A 


detailed description is given only for those reactions relevant to the production of guanine 
and xanthine. For a detailed description of all other reactions, see Michaelian [26]. 


1. 
2. 


10. 


11. 


12. 


Hydrogen cyanide HCN (H) hydrolyses to formamide HzNCOH (F) [5,91,94]; 

A photon-induced tautomerization coverts formamide (F) into formimidic acid 
(Fa) [55,56,90,92,93]; 

Formimidic acid (Fa) can, in turn, be photolysed into HCN (H) (or HNC) plus H20 
(dehydration) [92,93,97]; 

Formamide (F) hydrolyses to ammonium formate (Af) [93,94]; 

HCN (H) thermally polymerizes into (HCN),. Its most stable tetramer (x = 4) is 
known as cis-diaminomaleonitrile, cis-DAMN (C) [5]; 

HCN (H) can also thermally polymerize into trans-diaminomaleonitrile, trans-DAMN 
(T) [5]; 

Trans-DAMN (T) and cyanogen (Cg) are good catalysts for the polymerization of 
4HCN into cis-DAMN (see Table 6 of Sanchez et al. [5]). The catalytic effect of trans- 
DAMN on the tetramization of HCN was incorporated into the model by reducing the 
energy of the activation barrier such as to give the same amplification factor of 12 due 
to the catalytic effect of the inclusion of 0.01 M trans-DAMN in the HCN solution 
observed in the experiments of Sanchez et al. [5] at a temperature of 20 °C. 


The catalytic effect of cyanogen (Cg) is taken to be the same as that of trans-DAMN (T) 
since Sanchez et al. determined both of these to be strong catalysts [5]. As mentioned 
above, cyanogen is a precursor needed for guanine obtained from the Lyman-a line 
(121.6 nm) on atmospheric HCN giving, through photolysis, the CN~ radical with 
high quantum yield (reaction R177 Table 1 of Zahnle [51]), and this can form cyanogen 
(CN)2 by interacting with a second HCN molecule (reaction R199, Table 1 of Zahnle); 
Trans-DAMN also acts as an auto-catalyst for its own thermal production from 4HCN. 
Cyanogen is also an effective catalyst for this reaction (Table 6 of Sanchez et al. [5]); 
(a) By absorbing a photon of 298 nm, cis-DAMN (C) transforms (through rotation 
around a double carbon covalent bond) into trans-DAMN (T). (b) There is also a 
smaller quantum efficiency for converting trans-DAMN back into cis-DAMN by 
absorption of a photon at 313 nm [5,58,60]; 

A photon at 313 nm electronically excites trans-DAMN (T) which then transforms into 
2-amino-3-iminoacrylimidoyl cyanide, AIAC (J), through proton transfer from one of 
the amino groups [60]; 

AIAC (J), on absorbing a photon at 275 nm, transforms through photon-induced 
cyclicization (ring closure) into an azetene intermediate ((5) of Figure 3) in an electronic 
excited state, which then transforms to the N-heterocyclic carbene ((6) of Figure 3) and 
finally this tautomerizes to give the imidazole, 4-aminoimidazole-5-carbonitrle, AICN 
(I) ((7) of Figure 3) [60]; 

The imidazole AICN (I), created in the previous photochemical reaction #11, is con- 
verted through hydrolysis to 4-aminoimidazole-5-carboxamide, AICA (L) [6]; 
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13. Interaction of the imidazole AICN (I) with ammonium formate (Af), together with the 
catalytic effect of formamide (F), leads to adenine (A) [98,99]; 

14. Adenine (A) can also be obtained through the attachment of HCN (H) to AICN (I) to 
form amidine (Am), which is a formamide (F) catalyzed thermal reaction involving 
formimidic acid (Fa) [100]; 

15. A subsequent tautomerization of amidine (Am) is required (calculated to have a high 
barrier of about 50 kcal mol!) which, once overcome by absorbing a photon at 250 nm, 
allows the system to proceed through a subsequent barrier-less cyclicization to form 
adenine (A) [61]; 

16. Hydrolysis of adenine (A) gives hypoxathine (Hy), determined by Zheng and Meng 
to have a transition state barrier of 23.4 kcal mol~! [96]; 

17. The combination of AIAC (L) with cyanogen (Cg) (a precursor produced in the middle 
atmosphere from HCN—see reaction # 7) through a thermal reaction leads to guanine 
(G) and HCN (H) [6]; 

18. Hydrolysis of guanine (G) leads to xanthine (Xa) [89,95]; 

19. Photochemical reactions 19 to 28. These represent the absorption (within in a 20 nm 
region centered on the wavelength of peak absorption) and dissipation through inter- 
nal conversion at a conical intersection to the ground state on a sub-picosecond time 
scale. All molecules listed in this set of photo-reactions are photo-stable because of this 
peaked conical intersection connecting the electronic excited state to the ground state. 
These reactions, with large quantum efficiencies, represent the bulk of the flow of en- 
ergy from the incident UV-C spectrum to the emitted outgoing infrared ocean surface 
spectrum and therefore contribute most to photon dissipation, or entropy production. 


To simplify the kinetic equations for the photochemical reactions listed in Table 2, 
we assume that the molecules only absorb within a region +10 nm of their maximum 
absorption wavelength Ajay and that this absorption is at their maximum molar extinc- 
tion coefficient €may (Table 1), and finally that these wavelength regions do not overlap 
for calculating the shadowing effect of molecular concentrations at higher positions in 
the vesicle. 

We assume that the vesicle is at the ocean surface and the depth coordinate is divided 
into 1 = 20 bins of width Ax = 5 um and the time step for the recursion calculation for the 
concentrations is 10 ms. The recursion relation for the factor of light intensity L) (i,C) for a 
concentration C of the molecule, at a depth x(i) = i- Ax below the ocean surface will be, 


bCU)) = 14G=1,CG=1))e 10 oe, (1) 


where «, is the absorption coefficient of water at wavelength A and €) is the molar extinction 
coefficient of the particular absorbing substance which has concentration C(i) at x(i). 

The kinetic equations giving the increments in concentration after each time step 
At = dt, for use in a discrete recursion relation, at a depth x below the ocean surface are 
determined from the reactions listed in Table 2 to be the following: 
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where the differentials are calculated discretely (e.g., dH /dt = AH/At) and all concentra- 
tions are calculated at discrete time steps of At = 10 ms and the calculated value of the 
change (e.g., AH(j) /Af) for time step j is summed to the previous value (e.g., H(j — 1)). The 
day/night factor d is equal to 1 during the day and 0 at night. [920, In9g, 1313, 1275 and [50 are 
the intensities of the photon fluxes at 220, 298,313, 275 and 250 nm respectively (Figure 1). 
€, are the coefficients of molar extinction for the relevant molecule at the corresponding 
photon wavelengths 1. 


3.4. Vesicle Permeability and Internal Diffusion 


The permeability of the fatty acid vesicle will decrease both with increasing size of 
the molecule and with increasing size of its electric dipole moment, and increase with 
temperature. The final and intermediate product molecules have large sizes and large 
dipole moments compared with the precursors, implying a tendency towards entrapment 
within the vesicle. We assume that the vesicle is completely permeable to H2O, HCN (H) 
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formimidic acid (Fa) and cyanogen (Cg) but impermeable to all the other intermediate 
and final products. Note that ammonium formate (produced within the vesicle in reaction 
#4) would be in its ionic form and therefore also unable to cross the fatty acid membrane 
(permeability across lipid membranes are reduced by orders of magnitude if the molecules 
are polar or charged [101)). 

The diffusion constant Dy for the molecule Y within the inner aqueous region of the 
vesicle will depend on the interior solution viscosity, which is dependent on the concen- 
tration of organic material. Studies of inter-cellular diffusion of nucleotides indicate three 
factors influencing diffusion rates, besides temperature, at high solute densities; the viscos- 
ity of the solvent, collision rate dependent on solute concentration, the size of the molecules, 
and the binding interactions between molecules [102]. The diffusion constant of adenine in 
pure water has been determined to be D4 = 7.2 x 10-© cm? s~! [103] while the measured 
diffusion rates in the cytoplasm of a cell can be as low as 1.36 x 107° cm? s~! [102]. 

Films of organics and trace metals, with a high density of lipids and other hydrocar- 
bons, produced, for example, by the ultraviolet spectrum of Figure 1 on CO? saturated 
water [23], would be expected on the ocean surface during the Archean. Diffusion constants 
within this sea surface microlayer would then be significantly smaller than for bulk water. 
Diffusion rates inside the vesicle will depend on the amount of organic material captured at 
the air/water interface during the formation of the vesicle (which would have considerable 
spatial variability) and on the amount of ongoing organic synthesis within the vesicle. 

The Stokes-Einstein diffusion constant of uncharged spherical molecules is given by: 


ne kpT 
67tyr 


(17) 


where T is the absolute temperature, 7 the dynamic viscosity of the solvent, and r the radius 
of the particle (assumed spherical). For polar solvents such as water, molecules with a 
charge or dipole moment experience further drag due to the charge-dipole or dipole-dipole 
interaction [104]. All diffusion constants can, therefore, be approximated relative to that 
measured for adenine at 300 K through the formula: 


_ waar? T 


D 
a (2+ py) Al/? 300 


D,(300), (18) 


where Ay is the polar surface area of adenine, and 4 and ply are the dipole momenst of 
adenine and the molecule Y respectively (Table 1). The factor of 2 in the equation weighs 
the importance of including the dipole moment in the diffusion calculation (Equation (17)), 
which is complex but relatively unimportant to final results. 

Our analysis is performed with the diffusion constant for adenine found in present 
day cytoplasm [102] (i-e., D,(300) ~ 1 x 10~° cm? s~!). Using the values given in Table 1 
for the molecular dipole moments and their polar areas, and Equation (18), in Table 3 we 
list the results for the diffusion constants of all relevant molecules at 80 °C (353.15 K) with 
respect to that for adenine at 300 K, D, (300). 


Table 3. Diffusion constants at 80 °C (353.15 K) relative to that of adenine at 300 K (ie., 
dy = Dy(353.15)/D,(300) obtained from Equation (18). In this paper, D4(300) is taken to be 
1x 10° om’ s+. 


dc dr dy dy dy dam da dny dcg dg dxa 


3.84 2.26 


1.11 2.81 2.81 1.86 1.66 1.18 1.18 2.16 6.76 1.28 1.55 


Cyclical boundary conditions are assumed for diffusion, except for HCN (H), formimidic 
acid (Fa) and cyanogen (Cg) which, because of their small size and small dipole moment, 
can permeate the vesicle wall and, therefore, at the wall they are given their fixed value 
specified in the initial conditions of the environment outside the vesicle (see below). The 
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second order derivatives for calculating the diffusion were obtained using the second order 
finite difference method with double precision variables. 


3.5. Initial Conditions 


Miyakawa, Cleaves and Miller [91] estimated the steady state bulk ocean concentration 
of HCN at the origin of life assuming production, through electric discharge on a neutral 
atmosphere containing some methane, of radicals which attack N2, leading to an input rate 
to the oceans of 100 nmole cm~? y~!, and loss of HCN due to hydrolysis and destruction at 
submarine vents with a 10 million year cycling time of all ocean water for an ocean of 3 Km 
average depth. For an ocean of pH 6.5 and temperature of 80 °C, they obtained a value of 
[HCN] = 1.0 x 107!° M [91]. 

Not included in the Miyakawa et al. calculation, however, is, as mentioned above, 
the fact that HCN can also be produced through the solar Lyman alpha line (121.6 nm) 
photo-lysing N>2 in the upper atmosphere giving atomic nitrogen which then combines 
with CH and CH) to give HCN, or through 145 nm photolysis of CH, leading to a CH* 
radical which attacks N2 to give HCN [47]. Including this UV production would increase 
the input of HCN concentration to the oceans by a factor of at least 6 [51,105,106]. More 
importantly, the first ~100 um of the ocean surface is now known to be a unique region 
(the hydrodynamic boundary layer) where surface tension leads to enriched organics 
with densities up to 10+ or 10° times that of organic material in the water column slightly 
below [107]. Langmuir circulation, Eddy currents, and the scavenging action of bubbles, 
tend to concentrate organic material into this surface film. Trace metal enhancement in this 
microlayer can be one to three orders of magnitude greater than in the bulk [107,108]. 

This high density of organic material trapped through hydrophobic and ionic inter- 
actions at the ocean surface leads to significantly lower diffusion rates within the surface 
microlayer as compared to the ocean bulk [107]. Little diffusion and turbulence imply little 
mixing. The ocean microlayer is, therefore, a very stable layer which would not be recycled 
through ocean vents and thus allowed to accumulate. Finally, although HCN is very soluble 
in bulk water, molecular dynamic simulations have shown that it concentrates to about 
an order of magnitude larger at the air-water interface due to lateral HCN dipole-dipole 
interactions, and that HCN evaporates at lower rates than water [109]. 

Therefore, rather than assuming the low HCN bulk concentrations of Miyakawa et al., 
we instead consider a higher initial surface concentration of HCN (H) and cyanogen (Cg) 
of 6 x 10~° and for formamide (F) and formimidic acid (Fa) of 1 x 107°, the latter resulting 
from a photochemical tautomerization of formamide, the hydrolysis product of HCN 
(reactions #1 and #2 of Table 2). We also allow for perturbation of the system by considering 
the probable existence of small and sparse patches of much higher concentration of the 
precursors, justified by the characteristics of the ocean microlayer and the dipole-dipole 
interaction between HCN molecules. The system is perturbed by a single patch of rela- 
tively high concentration (0.1 M) of HCN, cyanogen (Cg), and formimidic acid (Fa) (see 
Michaelian [26] for a description) into which our vesicle is assumed to drift (at day 10.4) 
only once for 2 min during 60 Archean days. The initial concentrations of all other reactants 
inside the vesicle (assumed impermeable to these) are taken to be 1.0 x 10-19 M. 

Our system is under a diurnal 8 h constant flux of radiation followed by an 8 h period 
of darkness during which thermal reactions proceed, but not photochemical reactions. 
Different surface UV-C light intensities will be considered due to the uncertainty in the 
Archean atmospheric gas concentration profile. The “nominal” intensity corresponds to a 
weakly reducing and self-shielding (for A < 220 nm) atmosphere as proposed by Sagan 
and Chyba [30], which assumes that most of the greenhouse warming (required to solve 
the Faint Young Sun paradox [30]) is provided by CH, and a small amount of ammonia 
(mixing ratio 1 x 10~°) which is shielded by the products resulting from CO and CH, 
photolysis (e.g., HCN and aldehydes). The corresponding integrated (210-285 nm) UV-C 
light intensity for this “nominal” scenario (Fig. 1) is 4 W m~*. On the other hand, Cnossen 
et al. [32] suggest a neutral atmosphere in which the partial pressure of CO> has to be very 
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large, 0.24 bar, in order to achieve the required global warming with the faint young sun. 
This results in a much reduced integrated (210-285 nm) UV-C surface flux of 0.10 W m~?. 
We therefore perform our simulations for a number of different incident photon intensities 
covering both atmospheric profile scenarios. 


4. Results 
Evolution of the Concentration Profile 


In Michaelian [26] we have validated our photochemical reaction—diffusion model by 
comparing results to experimental data for the time dependence of a number of different 
intermediate products. 

In Figures 4—7 we present the concentrations of the relevant molecular products within 
the vesicle as a function of time in Archean days (16 h) in the dissipative photochemical 
synthesis of adenine and guanine, obtained by solving simultaneously the differential 
kinetic Equations (2)—(16), for the initial conditions and diffusion constants listed in the 
figure captions. 

The concentration profiles of the molecules evolve over time because of accumulation 
of photo-products within the vesicle and due to a deliberate external perturbation of the 
non-linear system at 10.4 days (0.1 M of HCN, cyanogen, and formimidic acid) which leads 
it to a new stationary state in which conversion of the environmental precursor molecules 
HCN and cyanogen into adenine and guanine occurs at greater rates. This leads to greater 
photon dissipative efficacy of the system, i.e., to a concentration profile of the product 
molecules which dissipates more efficiently the incident UV-C spectrum. 

In Figure 6 we plot the concentration profiles of the photochemical products obtained 
also at 80 °C but with 1/100 of the incident light flux (i.e., UV-C flux 1.25 W m~? m7! at 
260 nm, integrated flux 210-285 nm = 0.040 W m~*). This would correspond approximately 
to the illumination on Earth’s surface for a very neutral atmosphere if the only greenhouse 
gas to resolve the Faint Young Sun paradox (and increase the surface temperature to about 
80 °C) was CO, (see Section 3.2), assuming the standard solar model. It is interesting to 
note that the adenine concentration after 60 days actually increases slightly with this lower 
surface UV-C light flux, while all other purine concentrations decrease. This is because 
less UV-C light implies more formamide (F) being available (less of photoreaction #2) and 
thus more ammonium formate (Af) (reaction #4). Formamide (F) is a catalyst in the most 
important thermal reaction #13 which has a low activation energy barrier for converting 
AICN (I) and ammonium formate (Af) into adenine (A). 

Our resulting adenine concentrations are between about one and two orders of mag- 
nitude greater than those of guanine, which is consistent with experimental data of Levy 
et al. [10] using ammonium cyanide as a precursor under long time eutectic freezing con- 
ditions. This greater difficulty of producing guanine abiotically through both processes, 
eutectic freezing and photochemical dissipation, may explain why the major energy storage 
molecule of life is adenine triphosphate (ATP) rather than guanine triphosphate (GTP). 

Given the fixed concentrations of HCN (H), cyanogen (Cg) and formimidic acid (Fa) in 
the environment, to which the vesicle is permeable, photochemical reactions occur during 
daylight hours (marked as violet colored sections of the horizontal dashed line labeled 
as D/N). This gives rise to the observable diurnal oscillations in the concentrations of 
trans-DAMN (T) and AICN (I) since these are direct products of photochemical reactions. 
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Figure 4. Concentrations as a function of time in Archean days (16 h) of the precursor and product 
molecules; H—HCN, F—formamide, Fa—formimidic acid, Ce—cyanogen, Af—ammonium for- 
mate, C—cis-DAMN, T—trans-DAMN, A—adenine, I—AICN, J—AIAC, L—AICA, Am—amidine, 
Hy—hypoxanthine, G—guanine, Xa—xanthine, dissipatively structured on route to the synthesis of 
adenine (black trace) and guanine (blue trace). The initial conditions are; UV-C flux 125 W m-? um! 
at 260 nm (Figure 1), integrated flux 210-285 nm = 4 W m”, temperature T = 90 °C, initial con- 
centrations [H]p and [Cg]p = 6.0 x 10-5 M, [F]p and [Fa]y = 1.0 x 10-5 M and all other initial 
concentrations [Y]p = 1.0 x 107° M. The diffusion constant factor D 4 (300) was 1.0 x 10° cm? s~!. 
There is one perturbation of the system corresponding to the vesicle floating into a region of HCN 
(H), cyanogen (Cg), and formimidic acid (Fa) of concentration 0.1 M each for two minutes at 10.4 
Archean days. A new stationary state at much higher adenine and guanine concentrations is reached 

after the perturbation. The violet horizontal dashed line labeled “D/N” identifies alternate periods of 
daylight (violet) and night (blank). After 60 Archean days, the concentration of adenine within the 

vesicle (black trace) has grown by more than six orders of magnitude, from 1.0 x 10-1 to 1.19 x 10-4 

M and guanine (blue trace) by five orders of magnitude to 3.20 x 10-5 M. 
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Figure 5. The same as for Figure 4 except for a temperature of 80 °C. The adenine concentration 
reaches 2.4 x 1075 M and a guanine concentration of 7.6 x 10-7 M after 60 Archean days. 


At 10.4 Archean days, the vesicle is perturbed by assuming it passes through a region 
of high concentration of HCN (H), cyanogen, and formimidic acid (Fa), each of 0.1 M, for 
a 2 min period (vertical lines at 10.4 days). This sudden impulse in HCN, cyanogen, and 
formimidic acid concentration gives rise to rapid increases in all concentrations within the 
vesicle, in particular formamide (F), the hydrolysis product of HCN, which is an important 
catalyst for reaction #13 which produces adenine (A) from AICN (I) (see Table 2) and this 
reaction has a low activation barrier energy. Ammonium formate (Af) is used up in this 
reaction so its concentration decreases rapidly after the perturbation. More importantly, 
however, immediately after the perturbation there is a greater production of trans-DAMN 
(T) in the vesicle and since T acts as a catalyst for the polymerization of HCN (H) (reactions 
#7 and #8), this produces greater metabolism of HCN into DAMN within the vesicle, 
leading to a stronger diffusion of HCN into the vesicle from the outside environment. 
In other words, the reason that a short impulse of HCN, cyanogen, and formimidic acid 
give rise to an important increase in the rate of production of adenine and guanine is 
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that the vesicle’s semi-permeable wall, together with the set of equations describing the 
photochemical and chemical reactions, Equations (2)-(16), form a non-linear (auto- and 
cross-catalytic) system with multiple solutions at any given time. 
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Figure 6. The same as for Figure 5 (i.e., temperature of 80 °C) but with 1/100 of the incident light 
flux (i-e., UV-C flux 1.25 W m~? pm! at 260 nm, integrated flux 210-285 nm = 0.040 W m7“). This 
would correspond to the illumination on Earth’s surface if the only greenhouse gas to raise the 
surface temperature was CO2, assuming the standard solar model. The adenine concentration reaches 
3.8 x 10-5 Manda guanine concentration of 2.9 x 10-8 M after 60 Archean days. 


The perturbation at 10.4 days causes the system to leave the attraction basin of one 
solution determined by its initial conditions and evolve towards a different, and more 
probable, stationary state with a higher rate of production of adenine and guanine (slope of 
the black and blue traces respectively of Figure 5). The new stationary state is more stable 
than the initial state because its concentration profile is more dissipative, i-e., with more 
molecules having conical intersections dissipating the absorbed photon energy rapidly 
into heat, and, therefore, there is less probability of photochemical reactions reverting the 
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concentration profile to the initial state. There is thus a thermodynamic driving force for 
evolution to the new stationary state after the perturbation due to the system obtaining a 
greater photon dissipative capacity [26]. For systems where local thermodynamic equilib- 
rium (in space and time) can be defined (see Section 3 of Michaelian [26]), this implies that 
the entropy production of the system increases (see below). This example of the dissipative 
structuring a UV-C pigment concentration profile, elucidates the physics and chemistry 
behind a thermodynamic selection relevant to biological evolution at its earliest stages [26]. 
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Figure 7. The same as for Figure 4 except for a temperature of 70 °C. The adenine concentration 
reaches 4.3 x 10-© Manda guanine concentration of 7.5 x 10-8 M after 60 Archean days. 


The temperature dependence of the concentrations of the product molecules obtained 
after 60 Archean days is given in Figure 8. Ammonium formate (Af) is produced by 
the hydrolysis of first HCN (H) to formamide (F) (reaction #1) and then hydrolysis of 
formamide to Af (reaction #4). Both of these reactions have high activation energies and 
this results in Af only being produced in important quantities at temperatures greater 
than about 80 °C. Most of the adenine (A) production occurs through reaction #13 which 
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consumes Af and therefore high temperatures are important to the production of adenine. 
The hydrolysis of AICN (I) to AICA (L) (step (7) to (10) of Figure 3) also has a high activation 
energy barrier (reaction #12), also implying high temperatures necessary for the production 
of hypoxanthine (Hy), guanine (G) and xanthine (Xa), as can be seen by the slopes of their 
traces on Figure 8. 
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Figure 8. The temperature dependence of the concentrations of the product molecules af- 
ter 60 Archean days, with the initial conditions; UV-C flux 125 W m2 um! at 260 nm 
(Figure 1), integrated flux 210-285 nm = 4 W m7”, [H]p = 6x 1075, [Cg]lo = 6 x 10-%, 
[F]o = 1x 10-5, [Fa]o = 1 x 10-5 M, and all other molecules [Ylo = 1 x 10-19 and the diffusion 
constant D4 (300) = 1.0 x 10-©, with one perturbation of [H], [Cg] and [Fa] to 0.1 M for 2 min at 
10.4 Archean days. 


In Figure 9 we plot the concentrations of the product molecules obtained at a tem- 
perature of 80 °C after 60 Archean days as a function of the incident light intensity in 
the 210-285 nm region, by scaling the spectrum given in Figure 1. Maximum adenine 
concentration is obtained for fluxes of about 1.25 W m~? um! at 260 nm, however, at 
this flux, the concentration of guanine is 3 orders of magnitude smaller. At the nominal 
intensity of 125 W m~* um~! at 260 nm (Figure 1), the guanine concentration is about 
30 times lower than that of adenine which remains high at 2.4 x 107° M. 
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Figure 9. The dependence of the concentrations of the product molecules after 60 Archean days on the 
incident UV-C light intensity, with the initial conditions, T = 80 °C, [H]o = 6 x 10-5, [Cg]o = 6x 10-5, 
[Flo = 1x 10-5, [Fa]lo = 1x 10~°M, and all other molecules [Y]p = 1 x 10~!9 and the diffusion 
constant D (300) = 1.0 x 10~®, with one perturbation of [H], [Cg] and [Fa] to 0.1 M for 2 min at 
10.4 Archean days. 


In Figure 10 we plot the entropy production as a function of time in Archean days 
due solely to the absorption and dissipation of photons into heat through the conical in- 
tersections of the product molecules associated with reactions 19 to 28 of Table 2 for the 
particular concentration profile at the given time. For a description of the entropy pro- 
duction and its evolution through time in relation to the non-equilibrium thermodynamic 
Glansdorff-Prigogine universal evolutionary criterion, see reference [26]. 
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Figure 10. The entropy production as a function of Archean days for our UV-C photochemical dissi- 
pative structuring system within a vesicle floating at a sea surface at a temperature of 80 °C leading 
to the purines. Entropy production increases monotonically as photochemical reactions convert HCN 
into the photon dissipative product molecules, including adenine and guanine. During the day, 
entropy production is due to the dissipation of the UV-C spectrum into heat by the corresponding 
product concentration profile. During the night, entropy production is due to thermal chemical 
reactions, however, this entropy production is small and not included in the figure. At 10.4 Archean 
days, the system is perturbed and moves into a new stationary state with the entropy production 
increasing by more than 4 orders of magnitude and continues to increase as the concentration profile 
evolves (see Figure 5). 


5. Discussion 


As Figure 8 illustrates, high temperatures are very important to a building up high con- 
centrations of purines within the vesicle. High temperatures also foment phosphorylation 
with phosphate salts using formamide as a catalyst, favoring the formation of acyclonu- 
cleosides and the phosphorylation and trans-phosphorylation of nucleosides which only 
occurs efficiently at temperatures above 70 °C [83,84]. High temperatures are similarly 
important for the abiotic production of amino acids and fatty acids. Such high temperatures 
are, in fact, consistent with the geochemical data of the early Archean [70-74] and with the 
genome comparison of strongly conserved ribosomal RNA [75]. 

Important concentrations of all purines buildup within the vesicle within only about 
60 Archean days under the nominal environmental conditions considered here. There 
is no need to begin with unrealistically large initial concentrations of HCN by invoking 
eutectic freezing, and there is no need for highly alkaline conditions in order to favor HCN 
polymerization over hydrolysis. 
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It is instructive to compare our overall non-equilibrium results obtained with the 
model of UV-C dissipative structuring of adenine and guanine from HCN within a lipid 
vesicle at 80 °C with the quasi-equilibrium experiments of Ferris et al. [39] using high HCN 
concentration, alkaline pH, and room temperature. Starting with a concentration of HCN 
(0.1 M) in water (pH 9.2), and allowing this solution to polymerize in the dark at room 
temperature for a 7 months, and then subjecting these polymers to hydrolysis at 110 °C 
for 24 h, Ferris et al. obtain an adenine yield of 1 mg L7!, equivalent to a concentration 
of 7.4 x 10~° M (the molar mass of adenine being 135.13 g mol~!). Our model gives a 
similar adenine concentration after only 26 Archean days (Figure 5), starting from a much 
lower and more realistic initial concentration of HCN of only 6.0 x 10-5 M (with only one 
perturbation of 0.1 M for two minutes) and a more probable Archean neutral pH of 7.0 
at 80 °C and under a UV-C flux integrated from 210-285 nm of about 4 W m~? during 
daylight hours (Figure 5). At 90 °C our adenine concentration within the vesicle surpasses 
Ferris et al. value after the perturbation at only 10.4 days (Figure 4). 

As detailed in Michaelian [26], besides entrapment of product molecules inside the 
vesicle, another concentration mechanism for these molecules that may arise is the coupling 
between reaction and diffusion in the non-linear regime which leads to the breaking of 
spatial symmetry. For low diffusion rates, the homogeneous stationary state is no longer 
stable with respect to space dependent perturbations and intermediate and final products 
may become preferentially concentrated and be consumed in the center of the vesicle (see 
Michaelian [26]). 

The dependence of the final concentrations on the light intensity is complex (Figure 9) 
but important values of both adenine and guanine can be obtained for surface intensities cor- 
responding to the nominal value (for a weakly reducing atmosphere) of 125 W m-~* um7! 
at 260 nm (Figure 1) and for intensities either one order of magnitude smaller or greater. 
This result may have a dependence on the particular shape of the incident spectrum, but 
this was not studied here. 

Inorganic catalysts have not been included in our reaction scheme. These can increase 
significantly the rate of purine production. For example Cu*? ions increase the rate constant 
for the conversion of HCN (H) to cis-DAMN (C) [5] (reaction #5). Metal ions would have 
been in high abundance at the ocean surface microlayer [107,110]. 

The dissipative structuring of the pyrimidines under the same Archean UV spectrum 
(Figure 1) will be considered in a future article. 


6. Conclusions 


We have presented a simple model of UV-C photochemical dissipative structuring of 
the purines within a lipid vesicle permeable to the common precursor molecules H2O, HCN, 
cyanogen, and formimidic acid, but impermeable to the reaction products, floating on a hot 
Archean ocean surface. The chemical and photochemical kinetics are based on published 
experimental and quantum mechanical ab initio data. The surface physical conditions are 
consistent with biochemical and geochemical fossil evidence from the early Archean. 

High temperatures are important to obtaining high purine final concentrations which 
increase by roughly an order of magnitude for each 10 °C increase in surface temper- 
ature. Important concentrations of all purines buildup within the vesicle within only 
about 60 Archean days. There is no need to assume unrealistically high concentrations of 
HCN by invoking eutectic freezing or highly alkaline conditions. Our resulting adenine 
concentrations are between 1 and 2 orders of magnitude greater than those of guanine. 

The final product concentrations have a complex dependence on light intensity. The 
largest concentrations of adenine occur for light intensities of between 1 and 125 W m~? m=! 
at 260 nm, with guanine concentrations reaching about 1/30 of adenine concentration at 
125 W m~ um! at 260 nm (Figure 9). This adenine to guanine ratio is consistent with 
experimental data of Levy et al. [10] employing long time eutectic freezing conditions. 
This greater difficulty of producing guanine abiotically, through both eutectic freezing and 
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photochemical dissipation, may explain the fact that the major energy storage molecule of 
life is adenine triphosphate (ATP) rather than guanine triphosphate (GTP). 

Perturbations, caused by the vesicle floating into isolated patches of higher concen- 
tration precursors floating on the ocean surface, can provoke the non-linear auto-catalytic 
system into new states of higher product productivity with a greater “metabolism” of 
both precursor molecules (HCN and cyanogen) from the environment and the incident 
photons, leading to a concentration profile of much greater photon dissipative efficacy. 
This evolution towards stationary states of usually higher photon dissipative efficacy is the 
essence of biological evolution which continues today, but now involving the whole solar 
spectrum and the entire biosphere. 

Evolution during the Archean was towards concentration profiles of product molecules 
with an absorption maximum closer to the peak intensity of the incident UV-C spectrum 
at 260 nm (Figure 1) and with peaked conical intersections to internal conversion, both 
increasing the overall efficacy of the dissipation of the incident UV-C solar flux. 

For low diffusion rates, there can be significant coupling of reactions with diffu- 
sion, leading to non-homogeneous distributions of some of the intermediate products, 
in particular, with greater concentration of the products at the center of the vesicle (see 
Michaelian [26]). Such spatial symmetry breaking could facilitate further dissipative struc- 
turing, such as polymerization of the nucleobases into nucleic acid. 

Dissipative structuring, dissipative proliferation, and dissipative selection (as ex- 
plained in detail in Michaelian [26]) are the necessary and sufficient components for an 
explanation in physical-chemical terms of the synthesis, proliferation, and evolution of or- 
ganic molecules on planets, comets, asteroids, and interstellar space [18] and, in particular, 
for explaining the origin and evolution of the fundamental molecules of life on Earth. 
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Abbreviations 


The following abbreviations are used in this manuscript: 


AIAC 2-amino-3-iminoacrylimidoyl cyanide 
AICA 4-aminoimidazole-5-carboxamide 
AICN 4-aminoimidazole-5-carbonitrile 

CIT Classical Irreversible Thermodynamics 


DAMN _ diaminomaleonitrile 

DAFN diaminofumaronitrile 

PAHs Polyaromatic Hydrocarbons 

UV-A light in the region 360-400 nm 

UV-B light in the region 285-360 nm (only the region 310-360 nm is relevant here) 
UV-C light in the region 100-285 nm (only the region 210-285 nm is relevant here) 
UVTAR Ultraviolet and Temperature Assisted Replication 
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